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We consider Big Bang Nucleosynthesis (BBN) with long-lived charged massive particles. Before 
decaying, the long-lived charged particle recombines with a light element to form a bound state like 
a hydrogen atom. This effect modifies the nuclear reaction rates during the BBN epoch through 
the modifications of the Coulomb field and the kinematics of the captured light elements, which can 
change the light element abundances. It is possible for heavier nuclei abundances such as 7 Li and 
7 Be to decrease sizably, while the ratios Y p , D/H, and 3 He/H remain unchanged. This may solve 
the current discrepancy between the BBN prediction and the observed abundance of 7 Li. If future 
collider experiments find signals of a long-lived charged particle inside the detector, the information 
of its lifetime and decay properties could provide insights into not only the particle physics models 
but also the phenomena in the early Universe in turn. 

PACS numbers: 95.35.+d, ll.10.Kk, 12.60.-i 



INTRODUCTION 

Recent cosmological observations agree remarkably 
with standard ACDM models. The one and three- 
year data of Wilkinson Microwave Anisotropy Probe 
(WMAP) observation determined the cosmological pa- 
rameters to high precision [l], 0] . 

In light of such recent progress of cosmological obser- 
vations, it has been shown that the Universe should be 
close to flat, and most of the matter must be in the form 
of non-baryonic dark matter, which has been originally 
considered as one of the best candidates to explain an 
anomaly in the rotational curves of galaxies. 

In extension of the Standard Model explaining elec- 
troweak symmetry breaking and stability of the hier- 
archy, several candidates of the particle dark matter 
have been proposed such as the neutralino [3|, the grav- 
itino @, [EH, 0,11], the axino in supersymmetric the- 
ory, branon dark matter [Io| . Kaluza Klein dark mat- 
ter 11, 13] and Little Higgs dark matter 1^, 14\ and so 
on. The searches and the detailed studies of the dark 
matter have become one of the most exciting aspects of 
near future collider experiments and cosmological obser- 
vations. 

Considering such candidates in particle physics models, 
we expect the large amount of the dark matter particle 
will be produced at the near future colliders [15| . which 
will be powerful tools to understand the properties of the 
dark matter [16j . On the other hand, cosmological obser- 
vations may provide information in new particle physics 
models, and even some implications on undetectable the- 
oretical parameters in the collider experiments. Thus 
the connection of cosmology to collider physics may pro- 
vide wide possibilities to understand the properties of the 
dark-matter particle and check the cosmological models 



themselves. 

At the present stage, the detailed properties of the dark 
matter is still unknown. Therefore, even exotic proper- 
ties might be allowed. Future observations/experiments 
may prove them and single out or constrain dark-matter 
candidates. Even now, some problems in cosmological 
observations may already show some hints to understand 
the unknown properties of dark matter e.g., in the small 
scale structure problem 1?], HI l3| indicated in the cold 
dark matter halo, the low 7 Li problem [2C| and so on. 
There are several proposals to solve them by new physics 
M 0, 0, El EE US EH, HI M, However, 
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considerable astrophysical uncertainties may still exist. 

During the radiation dominated epoch well before the 
decoupling of the cosmic microwave background (CMB), 
it is not necessary that the dominant component of mat- 
ter is neutral, and that relic is the same as the present 
one. For stable CHArged Massive Particles (CHAMPs) 
33L 3{| , their fate in the universe had been discussed [34| , 
and the searches for CHAMPs inside the sea water were 
performed [36j , which obtained null results and got con- 
straints on stable CHAMPs [37]. According to their re- 
sults, the production of stable CHAMPs at future collider 
experiments is unlikely. However, such null results can be 
applied only for the stable CHAMPs, and still the win- 
dow for long-lived CHAMPs with a mass below O(TeV) is 
left open. Such possibilities for the long-lived CHAMPs 
were well-motivated in a scenario of super Weakly Inter- 
acting Massive Particle (super WIMP) dark matter 0], 
which may inherit the desired relic density through the 
long-lived CHAMP decays. The dominant component of 
the nonrelativistic (NR) matter during/after the BBN 
epoch might be charged particles. In supersymmetric 
theories, such a situation is naturally realized in grav- 
itino lightest supersymmetric Particle (LSP) and axino 
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LSP scenarios. Then the candidate for the long-lived 
CHAMP would be a charged scalar lepton [f| 0, HH . 

Trapping such long-lived CHAMPs, the detailed stud- 
ies of long-lived charged particle will be possible in future 
collider experiments, which may be able to provide some 
nontrivial tests of underlying theories, like measurement 
on the gravitino spin, on the gravitational coupling in 
the gravitino LSP scenario [391 ]. The trapping method 
in CERN LHC and International Linear Collider (ILC) 
has been performed in the context of supersymmetric 
theories [40( • Also the collider phenomenology [4lJ, |42j 
and the other possible phenomena 43j, |44| have been dis- 
cussed. 

In cosmological considerations of such long-lived par- 
ticles, the effects on BBN by the late-time energy in- 
jection due to their decays have been studied in de- 
tail [H, 0, E3, 0, Ei|. On the other hand, in the past 
studies of the effects on the light elements abundances, 
the analysis were simply applied to long-lived 'charged' 
massive particles, assuming all CHAMPs are ionized and 
freely propagating in the radiation dominated epoch well 
before the CMB decoupling. However, we show that 
these results are not always valid if the bound state with 
a CHAMP and light elements may have O(MeV) binding 
energy [351] . and the bound state might be stable against 
the destruction by the scattering off the huge amount 
of the background photons even during the BBN epoch. 
Also we show that heavier elements tend to be captured 
at earlier time. Namely the heavier light elements such as 
7 Li or 7 Be form their bound states earlier than the lighter 
light elements, D, T, 3 He and 4 He. Such a formation of 
the bound state with a heavy CHAMP may provide pos- 
sible changes of the nuclear-reaction rates and the thresh- 
old energy of the reactions and so on, which might result 
in the change of the light element abundances. 

What is the crucial difference from the case of electron 
captures? In case of the electron capture, since the Bohr 
radius of an electron is much larger than the typical pion- 
exchange length 0(l/m ff ), two nuclei feel the Coulomb 
barrier significantly before they get close to each other. 
On the other hand, in the case of the capture of the 
CHAMPs, the Bohr radius could be of the same order 
as the typical pion-exchange length. Then, the incident 
charged nuclei can penetrate the weakened Coulomb bar- 
rier, and the nuclear reaction occurs relatively rapidly. 
The importance of such a bound state in the nuclear re- 
action had been identified for cosmic muons [Hot 51 1. 
l 

Concerning a discrepancy in 7 Li between the standard 
big bang nucleosynthesis (SBBN) prediction by using the 
CMB baryon-to-photon ratio and the observational data, 



as we will show the details later, it is unlikely to attribute 
the discrepancy only to uncertainties in nuclear-reaction 
rates in SBBN P 



53l . l54j . However as we mentioned 



above, if CHAMPs exist, the nuclear reaction rates dur- 
ing the BBN epoch could be changed from the values 
known by experimental data or observations of the sun, 
and may potentially solve the current low 7 Li/H problem. 

If such long-lived CHAMPs existed and affected the 
light element abundances, the lifetime would be long (> 
lsec). They may be discovered as long-lived heavily ion- 
izing massive particles inside the detector in the collider 
experiments. The measurements of their lifetime and 
properties may provide new insights to understand not 
only the particle physics models but also the phenomena 
in the early Universe in turn. 

In this paper, we discuss the possible change due to the 
long-lived CHAMPs during/after BBN epoch and con- 
sider the effects on BBN. 2 



SBBN AND OBSERVED LIGHT ELEMENTS 

The theory in SBBN has only one theoretical param- 
eter, the baryon to photon ratio rj, to predict primor- 
dial light element abundances. Comparing the theoreti- 
cal predictions with observational data, we can infer the 
value of rj in SBBN. It is well known that this method 
had been the best evaluation to predict rj before WMAP 
reported their first-year data of the CMB anisotropy [l[ . 

WMAP observations have determined 77 in high preci- 
sion. The value of rj reported by the three-year WMAP 
observations Q is 



V 



(6.10 ±0.21) x 10 



-10 



(1) 



where rib is number density of baryon, and n 7 is number 
density of the cosmic background photon. In Fig [T] we 
plot the theoretical prediction of the light element abun- 
dances with their 2 a errors. The vertical band means 
the value of r\ reported by the three-year WMAP obser- 
vations at 2 a. 

We briefly discuss the current status of the theory 
of SBBN and the observational light element abun- 
dances below, and check the consistency with the CMB 
anisotropy observation. Further details of the observa- 
tional data are presented in a recent nice review by G. 
Steigman [55j |. The errors of the following observational 
values are at la level unless otherwise stated. Hereafter 
nx denotes the number density of a particle X . (X, C) 
denotes the bound state of CHAMP with an element X. 

The primordial abundance of D is inferred in the 
high redshift QSO absorption systems. Recently a new 



1 In muon catalysis fusion, the formation of an atom containing 
two nuclei may be important. 



2 In this paper, we use natural units for physical quantities. 



3 



data was obtained at redshift z — 2.525659 toward 
Q1243+3074 [H3|. Combined with these data [13, [H, [H, 
l60j . the primordial abundance is given as n.D/«-H|obs = 
(2.78+°;3g) x 10~ 5 . 3 It agrees excellently with the value 
of r] predicted in the CMB anisotropy observation. 

The abundance of 3 He can increase and decrease 
through the chemical evolution history. However, it is 
known that the fraction rt3 He /n,£> is a monotonically in- 
creasing function of the cosmic time [H, [6l[ . Therefore 
the presolar value is an upper bound on the primordial 
one, n m Jn D < 0.59 ± 0.54 (2a) 62}. In SBBN the the- 
oretical prediction satisfies this constraint. 

The primordial abundance of 4 He is obtained from the 
recombination lines from the low-metallicity extragalac- 
tic HII region. The mass fraction of the 4 He is inferred 
by taking the zero metallicity limit as O/H — > for the 
observational data [63j. A recent analysis by Fields and 
Olive obtained the following value by taking into account 
the effect of the Hel absorption, y(FO) obs = 0.238 ± 
(0.002) stat ± (0.005) syst , where the first and second errors 
are the statistical and systematic ones. On the other 
hand, Izotov and Thuan [64[ reported a slightly higher 
value, F(IT) obs = 0.242 ± (0.002) stat (±(0.005 WWiere 
we have added the systematic errors following [65ll66U67|. 
Olive and Skillman recently reanalyzed the Izotov-Thaun 
data [iH and obtained a much milder constraint [69j |. 
^(OS)obs = 0.249 ± 0.009. Even if we adopted the more 
restrictive value in Ref. [63|, SBBN is consistent with 
CMB. 

For 7 Li, it is widely believed that the primordial abun- 
dance is observed in Pop II old halo stars with temper- 
ature higher than ^6000K and with low metallicity as a 
"Spite'splateau" value. The measurements by Bonifacio 
et al. gave Logi [n7 Li /n H ]|obs = -9.66± (0.056) sta t± 
(0.06) sys . On the other hand, a significant dependence of 
7 Li on the Fe abundance in the low metallicity region was 
reported in [Til ■ If we take a serious attitude towards this 
trend, and assume that this comes from the cosmic-ray 
interaction [72| , the primordial value is 



|obs - (1.23_ 25 J x 10 

7lH 



(at 68% C.L). (2) 



Even if we adopt the higher value in Ref. [70(, the theo- 
retical prediction is excluded at 2 a outside the outskirts 
of observational and theoretical errors. Therefore when 
we adopt the lower value in ([2]), the discrepancy worsens. 
The centeral value of the observation is smaller than that 
of SBBN by a factor of about 3. This 7 Li problem has 
been pointed out by a lot of authors, e.g., see Ref. [20| . 

It has been thought optimistically that this discrep- 
ancy would be astrophysically resolved by some unknown 



systematic errors in the chemical evolution such as the 
uniform depletion in the convective zone in the stars . 
So far the researchers have added a large systematic er- 
rors into the observational constraint by hand [74], [75| . 

However, recently the plateau structure of 6 Li in nine 
out of 24 Pop II old halo stars was reported by Asplund et 
al. [7(| . The observed values of the isotope ratio na Li / m Li 
uniformly scatter between ~ 0.01 and 0.09 at 2 a, inde- 
pendently of the metallicity, and are approximately sim- 
ilar to the previous observational data (= 0.05 ± 0.02 at 
2 a jz3). Because the estimated 7 Li abundance in such 
stars is n7 Li /nH |obs = (1.1 — 1.5) x 10~ 10 , the upper bound 
on the primordial 6 Li agrees with SBBN. Although so far 
some models of the 6 Li and 7 Li production through the 
cosmic-ray spallation of CNO and a-a inelastic scatter- 
ing have been studied, the predicted value of ne^-Jm-n or 
nej, \ /n yi is obviously an increasing function of a metallic- 
ity El m> IH El . 

As we have discussed, to be consistent with the SBBN 
prediction and WMAP observations, we need a certain 
uniform depletion mechanism of 7 Li. Because 6 Li is more 
fragile than 7 Li, whenever 7 Li is destroyed in a star, 6 Li 
suffers from the depletion, too. If we require the pri- 
mordial abundance of 7 Li to be uniformly depleted to a 
smaller value by a factor of three, the ratio 6 Li/ 7 Li might 
have to be reduced by a factor of O(10) [12]. Therefore, 
we do not have any successful chemical evolution mod- 
els at the present, to consistently explain the observa- 
tional value of 6 Li/ 7 Li by starting from the theoretical 
prediction of the primordial values of 6 Li and 7 Li in the 
framework of SBBN. 

Thus, by adopting the r\ predicted in the CMB ob- 
servations, we would now have to check SBBN itself or 
modified scenarios related with BBN compared with the 
observational light clement abundances. 

In recent studies, it has been pointed out that the un- 
certainties on nuclear-reaction rates in SBBN never solve 
the discrepancy of 7 Li between the theory and the obser- 
vation. That is because the uncertainties are highly con- 
strained by known experimental data and observations of 
the standard solar model. In Ref. [H3|, the possible nu- 
clear uncertainties were investigated. It was shown that 
only a nuclear reaction rate more than 100 times larger 
in 7 Be(n,a) 4 He and 7 Be(d,p)2 4 He might provide sizeable 
change in the 7 Li abundance. Notice that 7 Be(n,a) 4 He 
does not have an s-wave resonance due to the symmetry 
of the outgoing channel while 7 Be(n,p) 7 Li has it. Since 
in the important energy region in the SBBN reaction 
T ~ 50 keV, which is near threshold of the processes, 
the contribution to 7 Be from 7 Be(n,a) 4 He is negligible 



3 Some of the observed data have large dispersion than ex- 4 See the recent report about spectroscopic observations of stars in 
pected and might have systematic errors which may cause higher the metalpoor globular cluster NGC 6397 that revealed trends of 

D/H 56, 5§| . atmospheric abundance with evolutionary stage of Lithium |73l 
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FIG. 1: Theoretical predictions of Y v , D/H, 3 He/H, 7 Li/H 
6 Li/H, 3 He/D and 6 Li/ 7 Li as a function of the baryon-to- 
photon ratio r\ in standard BBN with their theoretical errors 
at 95 % C.L. The WMAP value of rj at 95 % C.L.is also indi- 
cated as a vertical band. In the comparison between the BBN 
prediction and the central value of the observed abundances, 
it has been pointed out that the SBBN prediction with the 
WMAP value of the r\ shows too high by a factor of a few in 
7 Li abundance and too low by several orders of magnitude in 
6 Li abundance if there is no late time 6 Li production other 
than BBN. [3 



in SBBN relative to 7 Be(n,p) 7 Li because of the p-wave 
nature of the process. For 7 Be(d,p)2 4 He, the possibil- 
ity may not work in the light of the recent experimental 
Also Cyburt et.al. 



data 84]. Also Cyburt et.al. [83| discussed the uncer- 
tainties on the normalization of cross section for the pro- 
cess 3 He(cv,7) 7 Be and found that the uncertainties are 
constrained in the light of a good agreement between the 
standard solar model and solar neutrino data. 

Therefore the remaining possibilities may be uncer- 
tainties on the chemical evolution of Li from the BBN 
epoch to the present or effects due to new physics. Be- 
cause now we do not have any successful chemical evolu- 
tion models, it must be important to consider the effect 
of new physics. 

As we mentioned before, the existence of CHAMPs 
might provide possible change of nuclear-reaction rates 
during the BBN epoch, which may have some impact on 
the prediction of primordial light element abundances. 
In the next section, we will discuss the properties of the 
bound state and the recombination of CHAMP and the 
possible change of nuclear reaction rates. 



BOUND STATE WITH A CHAMP AND A LIGHT 
ELEMENT 

Evaluation of binding energy 

We evaluate the binding energy for the bound state of 
a negatively charged massive particle and a light element. 
We simply consider the case that the charged particle is 
a scalar. The extension to a fermion or the other higher 
spin cases would be straightforward although there exist 
little differences. Here we follow the way to evaluate 
the binding energy assuming uniform charge distribution 
inside the light element according to Ref. [35j . Then the 
Hamiltonian is represented by 



H = 



Z x Z c a 
2r x 



2m x 

for short distances r < rx, and 

Z x Zca 



Z x Z c a ( r 
2r x [ r x 



H 



V 
2m x 



(3) 



(4) 



for long distances r > r Xl where a is the fine struc- 
ture constant, r x ~ 1.2 A 1 / 3 / 200 MeV -1 is the nu- 
clear radius, Z x is the electric charge of the light el- 
ement, and Zc is the electric charge of the negatively 
charged massive particle. A is the atomic number, and 
m x is the mass of the light element X . Here we assumed 
m x <C mc ~ O(100 GeV), which means the reduced 
mass 1/fi = l/mc + I/tox ~ 1/mx- 

For large nuclei, the exotic charged particle may be 
inside the nuclear radius. The binding energy may be 
estimated under the harmonic oscillator approximation 
by 

3 ZxZca 1 Z x Z c a 

For small nuclei, the binding energy may be estimated 
well as a Coulomb bound state like a hydrogen atom, 



1 



bin 



-Z x z Zc 2 a 2 mx 



(6) 



For intermediate regions in between the above cases, by 
using a trial wave function, we can express 

Eu- m — — ( — - — F(ZxZcam x rxj), (7) 
UIli rx mxrx 

where F{x) is variationally determined 35]. For < 
ZxZcdmxTx < 1) the Coulomb model gives a good 
approximation. On the other hand, the harmonic os- 
cillator approximation gives a better approximation for 
2 < ZxZcctmxrx < oo. 

The binding energies are shown in Table U For a 
CHAMP with Zc = 1 and lighter elements (p, D, and 
T), typically ZxZcamxrx < 1- Thus the Coulomb ap- 
proximation works well. However, for heavier elements 
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Nucleus(X) 


binding energy (MeV) 


atomic number 


P 


0.025 


Z=l 


D 


0.050 


Z=l 


T 


0.075 


Z=l 


3 He 


0.270 


Z=2 


4 He 


0.311 


Z=2 


5 He 


0.431 


Z=2 


5 Li 


0.842 


Z=3 


6 Li 


0.914 


Z=3 


7 Li 


0.952 


Z=3 


7 Be 


1.490 


Z=4 


8 Be 


1.550 


Z=4 


io B 


2.210 


Z=5 



The thermal-averaged cross section is written as 



TABLE I: Table of the binding energies for the various nuclei 
in the case of Zc = 1 given in Ref . J35| • For elements heavier 
than 8 Be, the binding energies are given by the harmonic 
oscillator approximation. 



such as Li or Be, there may exist deviations which are 
more than O(10) percent. For elements lighter than 8 B, 
the binding energy is still below the threshold energy of 
any nuclear reactions. If the atomic number is not large 
like Li and Be, we can ignore the effects due to finite 
size and the internal structure (excitations to higher lev- 
els and so on) as a good approximation to calculate the 
capture cross section and the nuclear-reaction rates. 



CAPTURE OF CHAMPS IN THE EARLY 
UNIVERSE 

Recombination cross section 

We evaluate the recombination cross section from the 
free state to the IS bound state assuming a hydrogen- 
type bound state through a dipole photon emission [85| 
and a pointlike particle for the captured light element. 
Then the cross section is 



2 9 ir 2 aZ 2 x E hh 



■(- 



3 m x v K E hin + \m x v 2 ' 



2»w'aZ x E hir , 



3e 4 



(8) 



where v is the relative velocity of a CHAMP and a light 
element. Note that we have mxv 2 /2 ~ 3T/2 <C E^ in 
for NR particles in kinetic equilibrium. Here we use 
the Coulomb model (hydrogen type) to evaluate the 
capture rate [33|], where the binding energy Ebm = 
a 2 Z^,Z x mx/2 and the Bohr radius r^ 1 ~ aZcZxmx- 



{a r v 



-(ttttt) / d 3 Pl d 3 p 2 e- t (T r v 



nin 2 (2ir) 



' (ttAt) 2 /^" 



riQn r (2tt) 3 
>9_, 7 2 



2 J naZ x V2n E hin 

3^4 m 2 x y/mxT' 



where mc = mi+777,2 and u = mxmc / (mx+mc) 
with 



g t I d 3 p r a r ve 2 » T 



(9) 



m x 



n G 



(2tt) 3 
9 



(27T) 



d 3 pce T e 



d 3 p r e 



Here we have assumed that only one CHAMP is cap- 
tured by a nucleus. Since the photon emission from a 
CHAMP is suppressed, the recombination cross section 
for the further capture of an additional CHAMP by the 
bound state would be much smaller. Therefore, as a first 
step, it would be reasonable to ignore the multiple cap- 
ture of CHAMPs by a nucleus. 

Here we have estimated only the direct transition from 
the free state into the IS bound state. However, if the 
transition from higher levels into the IS state is suffi- 
ciently rapid against the destruction due to scatterings off 
the thermal photons, even the capture into the higher lev- 
els might contribute to the recombination of a CHAMP. 
The typical time scale of the transition from nth level into 
the IS state is l/{E h - m:ls - E hin:n ) ~ 0(1/E hin:13 ) where 
Ehuv.n is the binding energy of the nth level. Up to some 
levels, this time scale might be shorter than the destruc- 
tion rate after the IS state became stable. However, such 
higher-level captures would not significantly enhance the 
recombination cross section because the capture rate into 
higher levels is relatively suppressed and small. 

For highly charged massive nuclei or elements heav- 
ierthan boron, the binding energies with CHAMPs can 
become of the order of magnitude of the excitation ener- 
gies of nucleons inside the nuclei, or even of the same or- 
der of magnitude of the nuclear binding energies. In such 
cases, the capture process of light elements by CHAMPs 
may be nontrivial. In addition, to correctly calculate the 
capture rates, we would have to understand the modifi- 
cation by the effects due to not only the finite size but 
also the internal structure of the light element. In this 
paper, we ignore these effects because they are unimpor- 
tant since we consider lighter nuclei up to Li and Be. 



Case in kinetic and chemical equilibrium 

To evaluate the number density of the captured 
CHAMPs, we would be able to use the thermal relation 
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among chemical potentials if the capture reactions well 
establish the chemical equilibrium between the CHAMPs 
and the light elements. The number density is deter- 
mined by the following Saha equation, 



2 ,1ttT 3/2 



n^ 



rax 



(10) 



where nx and n 7 are number densities of a light element 
A and thermal photons, and E^i n is the binding energy 
of the light element. 



General cases 

However the question of whether such kinetic and 
chemical equilibrium are well established among all light 
elements and CHAMP is nontrivial. Here we consider 
the Boltzmann equations for CHAMPs, a light element 
A and the bound state (X,C). For CHAMPs, 



d_ 

dt 



n c + 3Hn c = 



d_ 

dt 



Re- 



capture 



(11) 



where H is the Hubble expansion rate. For a light ele- 
ment X, 



d_ 

dt 



n x + 3Hn x = 



d_ 
dt 



fusion 



+ 



d_ 

dt 



n x 



capture 



(12) 



For the bound state, 
d 



' d 




r d 


dt n ^ x \ 


fusion 


di nx _ 



capture 



(13) 



By using the detailed balance relation between the for- 
ward process X + C — > 7+ (X, C) and the reverse process 
(A, C) + 7 — > A + C, the capture reaction may be written 

by 



r d ' 




' d 




di nx _ 


capture 


dt nC _ 


capture 



■(a r v) [n c n x - n (C \X)n^(E > E hin )] , (14) 



where 



and 



n^E > E hin ) = n~ 



2C(3) v 2ttT 



2C(3), 



(15) 



(16) 



For a light element, if (a r v)n c /H 3> 1 is satisfied and 
the kinetic equilibrium is well established, we can get 
the Saha equation by requiring an equilibrium condi- 
tion [Jjfixlcapture = in this equation. Since we are 
interested in the time evolution of not only CHAMPs 
but also light elements, we carefully study the case of 
(a r v)nc / H > 1 even in the case of (<r r v)nx / H <C 1. 



Critical temperature at which a bound state is 
formed 



When the temperature is higher than the binding en- 
ergy of light elements, the destruction rate of bound 
states by scatterings off the thermal photons with E > 
Ebi n is rapid. Then only a small fraction of bound states 
can be formed, n^ c ,x) ~ n c n x /n 7 (E > E b i n ) <C n X - 
Once the temperature becomes lower than the binding 
energy, the capture starts, and the bound state becomes 
stable if the other destruction processes among the nuclei 
are inefficient. 5 The critical temperature at which the 
capture becomes efficient is estimated as follows. In the 
case of nx > nc, taking nc ~ ntc,x), we get a relation, 



(^f/V^^.OflO- 10 ). 



(17) 



On the other hand, in the case of nx < nc, taking nx 
n (c.x), we have 



('»*)3/2 c -^ 



Up 

n 



TOOGeV. . il c 



' O(10- 10 )(^^)(^-). (18) 
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m c 0.23' 



This analysis shows that the critical temperature is ap- 
proximately 



Ebin 

~40 



(19) 



In case of Z c = 1, we find T c - £ bi „/40 - 8kcV for 4 Hc. 

Here we consider the temperature where some frac- 
tion of A is captured by CHAMPs. For example, taking 
n{c,x)/n x ^ 1CT 5 , we get 



( m* )3 / 2e 



n c n x 



(2) 

This condition is satisfied at Tc ~ i?bin/30. Since the 
abundance of 4 He is large below O.lMeV, even though 
the only small fraction of 4 He is trapped by CHAMPs, 
there might be relevant effects caused by the captures. 

For protons, the efficient captures start at a temper- 
ature lower than 1 keV (at cosmic time longer than 10 6 
sec) . Since the bound state is neutral for a single-charged 
CHAMPs Zc = 1, and might be negatively charged for 
a multi-charged CHAMPs Zc > 1, there is no Coulomb 
repulsion anymore. Thus, even the bound state can col- 
lide with each other. If the number density of CHAMPs 



0(1Q- 6 ). (20) 



5 Note that the abundances of heavier elements such as Li and Be 
are smaller than those of lighter elements (p, D, T and He). As we 
will see later, considering the relic density of relevant candidates 
of CHAMPs, their capture can only affect on the abundance 
of the heavier elements. Our scenarios would not significantly 
change the lighter element abundances. 
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is not too small, and most CHAMPs are captured by 
protons, the change could be sizable for longer lived of 
CHAMPs (t > 10 6 sec). 6 



Capture rate 

Since the capture process competes with the expansion 
of the universe, we have to check if the following relation 
holds during the meaningful time, which ensures that the 
capture by CHAMPs is efficient compared to the expan- 
sion rate of the universe, 



H <C (<j r v)nc- 



(21) 



That is, the capture rate of a light element is controlled 
by the following k, 



(a r v)n c 




k is approximately 2.6 and 0.43 for 7 Li and 4 He at their 
critical temperatures, respectively. Here we assumed that 
n c ^ 0.23 and m c = lOOGeV. 

In the evaluation of the capture rates for light elements, 
we considered relatively large number densities of the 
CHAMPs, which are approximately similar to that of 4 He 
or even more because here we assumed that a CHAMP 
can decay into much lighter dark matter or almost mass- 
less SM particles later. Under these circumstances, we 
naturally expect a larger value of the capture rates than 
the upper limit in case of the stable CHAMP scenario. Of 
course, we have to check that the decay never disturbs 
the successful concordance of cold dark matter (CDM) 
with large scale structure formation in the Universe and 
so on. Later we will discuss this problem. 

Next let us estimate the time evolution of X itself and 
the capture fraction of A by a CHAMP. At below the 
critical temperature T c , the destruction term of {X, C) 
becomes negligible due to the Boltzmann suppression. 7 
Then the number densities of the light element X and the 
bound state of X with a CHAMP, (X, C) are obtained 
by solving the following equations. Here any destruction 
reactions of X would be negligible close to the end of the 
BBN epoch (< 50 keV.), 



-( 



((jy)n c r]x 
HT r, x (T c y 



(23) 



T/T 



FIG. 2: V(c,x)/Vx as a function of T/T c for 
(o r v)nc I H\t=t c =1, 3, 5, 10, and 30 from left to right, 
respectively. Here we have ignored the standard BBN 
processes. Also we have taken the initial condition as 
f?x(?c) ~ 0. If {a r v)nc /H\t-t c S> 1, the Saha equation will 
be a good approximation and the capture will immidiately 
occur at T ~ T c . On the other hand, if {a r v)nc / ' H\t=t c ~ 1 
or less than 1, the approximation by the Saha equation may 
be falied. 



d . V(c,x) > _ {<J r v)n c fix 
dT { r, x {T c )>~ 



HT r!x{T c ) 



(fusion part), 



where r)i — rii/s, and r]x{T c ) is the initial number den- 
sity per entropy density when the capture starts, assum- 
ing that the standard processes of the light elements are 
(almost) frozen out. We also assumed nc ^> nx which 
is correct except for 4 He. We find that if k is larger than 
unity at the critical temperature, the capture will be ef- 
ficient. 

Ignoring the fusion part of the standard processes 
in Eq . ([23]) . we find the following analytical solution of 
Vx(T), 



r lx (T) = r 1 x(T c )e- 2 ^ 1 -V^), 



(24) 



where k% = (a r v)nc / H\t=t - 8 For the numerical solu- 
tion of Eq. (JlU), see Fig. H 

For a more precise analysis, especially for the Boltz- 
mann equation of the CHAMP bound state, we may 
have to take into account the nuclear reaction processes 
simultaneously. Therefore we will need to do the numeri- 
cal calculations to solve the Boltzmann equations includ- 
ing both the capture and the BBN processes in the fu- 
ture [jjj ]. However, to qualitatively understand how large 



6 Since the CHAMPs with a long lifetime more than 2> 10 e sec 
may induce the other effects on cosmology [24| . 

7 The ignorance of the destruction term at T c may be valid if the 
recommbination cross section is not too large. If the cross section 
is enough large, the number density of bound state may be well 
descriobed by the Saha equation. 



For 7 Be, 7 Li and lighter elements, above approximation works 
well if Z c is close to 1. As we will see later, the change of nuclear 
reaction rates does not modify the fusion part of the noncapturcd 
light elements so much because the most of reverse processes 
has already been decoupled even after the other elements are 
captured by CHAMPs. 
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changes would be possible, for simplicity we assume only 
the instantaneous captures in the current work. 



CHANGE OF NUCLEAR REACTION RATES IN 
BBN BY THE CAPTURE OF CHAMPS 

The capture of light elements by CHAMPs weakens the 
Coulomb barrier in the nuclear reactions during/after the 
BBN epoch. The change of nuclear reaction rates could 
become large because the Coulomb factor exponentially 
suppresses the reaction rates. In general, the reaction 
rates among charged nuclei during the BBN epoch are 
determined by the competition between the coulomb sup- 
pression and the Boltzmann suppression, which play im- 
portant roles to determine the freeze-out of light element 
abundances at the end of the BBN epoch. Considering 
the corrections on these two exponential suppressions, we 
will next consider the possible changes of nuclear reaction 
rates. 



Coulomb potential and scattering problem 

If there are Coulomb expulsion forces, the wave func- 
tion of an incident particle would be exponentially sup- 
pressed at the target. Since we use a plane wave for 
the wave function to evaluate the incident flux at a suf- 
ficiently far place from the target, the real flux which 
is associated with the reaction would be evaluated by 
renormalizing the wave function. Since the change of the 
wave-function normalization from the plane wave is as- 
sociated with the state before the nuclear reaction, it is 
independent of the short distance nuclear reaction by nu- 
clei. We can expect that the Coulomb factor is factorized 
as follows 9 , 

2irZ a Z b a/v 2irZ a Z b a _ ^z a z ia 

Fab{v) = ^ ZaZha e « . (25) 

e v — 1 v 

After a CHAMP is trapped by a light element a, for 
a collision between a bound state (CHAMP + the light 
element a) and a light element 6, 



F {aC )b(P) 



2-KoiZ( aC) Z h 
' 



(26) 



9 This factorization may be valid only if the Bohr radius of bound 
state is not too large relative to the radius of the bound and 
incident nuclei. If the Bohr radius is large, which may be ex- 
pected in Zx = 1 nuclei cases, we have to understand how the 
bound state is disturbed by the incident nucleus. In such large 
Bohr radius cases, for example, to proceed nuclear fusion, the 
hydrogen-type bound state of the nucleus and a CHAMP may 
have to constitute a molecule before the nuclear fusion. Then we 
will have to evaluate the capture reaction rate of molecule. 



where z ( a C) = Z a — Zc- Note that (3 is the relative ve- 
locity between the bound state (aC) and the b element, 
not the a and the b element. Hence (3 could be slightly 
different from v which is the normal relative velocity be- 
tween the thermal a and the thermal b. Here we assumed 
that a light element can capture only one CHAMP (with 
the charge Zc)- 

For the case of nuclear reactions through a collision 
between charged bound sates (CHAMP + light ele- 
ment a and CHAMP + light element b), the Coulomb- 
penetration ability is determined by the relative velocity 
between the bound states. That is, 



(aC)(bC) 



—K , (27) 



where (i 2 = Pr / V( a c)(bC) - 0(T/m c ) « ft = 
0(T/mx)- Under these circumstances, the collision 
between charged bound states may be highly sup- 
pressed relative to the standard BBN reactions because 
z {aC) z {bC)/fc > Z a Z b /v if Z x > Z c . This bound state 
- bound state collision might become important if a huge 
number of CHAMPs are captured by 4 He. However, the 
typical temperature to start capture is below O(10) keV, 
and the Coulomb factors for the normal nuclear reactions 
in SBBN is highly suppressed, and have already been de- 
coupled by that time. Thus this type of collision will not 
contribute to any sizable changes of the light element 
abundances 

For Zx = 1 (Zc = 1) cases like protons, since there 
is no Coulomb suppression because the bound state is 
neutral, the collision between two bound states may be 
important. 



SBBN and thermal-averaged fusion rates 

First, we discuss nuclear-reaction rates in SBBN, and 
next we will extend the discussions to the cases with the 
CHAMPs. 

For simplicity, we consider the case of 2 — > 2 non- 
resonant reactions among charged light nuclei. The other 
cases may be straightforward through similar discussions. 
In a SBBN process a + b — > c + d, the forward process 
and the reverse process are defined by the difference be- 
tween the total masses in the initial and the final state. 
If Qab,cd = m a + m b -m c -m d = Qsbbn > 0, the process 
a + b — > c + d has no threshold and is called the forward 
process. On the other hand, the process c + d — > a + b has 
threshold (Q- value Q c d,ab = — <3sbbn < 0) and is called 
the reverse process of a + b ^ c + d . Usually the reverse 
process c+d — > a + b has a strong Boltzmann suppression 
by e -QsBBN/T Q value is larger than the Gamow 

peak energy of the process. 
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SBBN reaction rates with no threshold 

Naively the nuclear reactions of SBBN occur at almost 
the threshold region. Thus the cross section may be well 
described by the lower partial wave modes. Taking into 
account for the discussion of the wave function normal- 
ization in previous section, the reaction cross section is 
written as follow. 



(as + <T P v 2 + ....)F ab (v) 

, 2tt Z a Z b a ■2*z a z b 
a v(v) e 



(28) 



where oqv(v) = as + <Jpv 2 + 

Here we introduce a new variable, the "astrophysical S 
factor" which astrophysicists have used in the calculation 
of nucleosynthesis, 



S(E r ) — Ufusion-Erev e t 

= crav(v)TrZ a Z b a^ab 



(29) 



where E G = 2tt 2 Z 2 Z 2 a 2 fi ab and E r = pj/2fi ab = 
fiabV 2 /2. Notice that this S factor is a function of the 
center-of-mass (CM) energy and is inferred by the mea- 
surements of (Tfusion'f i n experiments and observations. 
The recent fitting functions are given in Refs. 5^, 53j |. 

By using this S factor, we calculate the thermal- 
averaged cross section. 




(x-x ) 2 + ....) 



3_ 1/3 

g 4I/3 G 



(30) 



where x = E r /T, x G = E G /T, and x = (xg/4) 1 / 3 . 
Since the main contribution of this integral comes from 
the stationary point of the exponent, we expanded the 
exponent around the stationary point xq — (xq/4) 1 / 3 . 

Finally we can evaluate the thermal-averaged nuclear- 
reaction rate among charged light elements. 



(Cfu 



32 E G /3 S(x T) 



4i/3 3 Mafc T 2 /3 
where l/fi a b = l/m a + l/m b . 



-e 4 



I E G U/3 

T > j ( 31 ) 



detailed balance relation. For example, in a 2 
resonant reaction a + b — > c + d, 



2 non- 



(q-fusion^)cd _ ^afe^3/2^ TO a + m b ^3/2 9a9b 

(o-fusionw)afc fJ-cd m c + m d g c g d 

where Q is the Q value of the forward reaction and g a is 
the number of degrees of freedom of the light element a. 
Notice that the factor e~®/ T arises from the Boltzmann 
suppression for the high energy component with E r > Q 
in thermal distribution. 



Extension to BBN with the captured CHAMP 

We have shown that the collision among charged 
CHAMP bound states will not result in any changes to 
SBBN. Here we focus on the nuclear- reaction rate for the 
collision between a bound state (CHAMP+light clement) 
and an unbound light element. 10 

Forward and Backward process 

Here we discuss the modifications of the short-distance 
nuclear- reaction rates mainly governed by the strong in- 
teraction. In CHAMP BBN (CBBN), the correspond- 
ing dominant process for the SBBN forward process 
a + b — > c + d may be (a, C) + b — > (c, C) + d or c + d + C, 
assuming (6, C) does not have sufficiently a large bind- 
ing energy against scattering of background photons, i.e., 
Ebin/T <C 40. Here (c, C) has a larger binding energy 
than that of (d, C). If the following condition is satisfied, 



QsBBN — ^bin aC > 0, 



(33) 



the final state is given by (a, C) + b — > c + d + C. 

On the other hand, even if the above condition is not 
satisfied, but if the following condition is satisfied, 



QcBBN — QsBBN + E] 



bin, cC 



Ebin,aC > 0, 



(34) 



(a, C) + b — > (c, C) + d is kinematically allowed, and the 
CHAMP in the final state will be trapped again. However 
if the bound state (c, C) does not have enough binding 
energy against the destruction due to thermal photons, 
the (c, C) state will be destroyed soon after the process, 
and the element c and the CHAMP will become free. 

For the Z G = 1 case and the relevant nuclei, because 
most of the Qsbbn values are sufficiently large, the case 
that Qsbbn > but (Jcbbn < would be rare. How- 
ever, in general, it might be possible. In such cases, even 



SBBN reaction rates with threshold 



We often evaluate reverse reaction rates from the ex- 
perimental data of forward reaction rates by using the 



10 For the case of scatterings among neutral bound states, the col- 
lision can easily occur. In such cases the calculation is straight- 
forward. 
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though the SBBN process does not have any threshold, 
the CBBN can have it. But the sign flip in the Q value 
occurs when the binding energy of a bound state with a 
CHAMP exceeds the nuclear binding energy of the pro- 
cess, which may mean that the bound CHAMP is not a 
spectator in the nuclear-reaction any more. In our fol- 
lowing analysis, we do not consider this kind of special 
cases. 

Next we simply assume that QsbbnOcbbn > 0. Let 
us consider the reverse processes ofa + 6 — ► c + din 
CBBN, which has a threshold characterized by <3cbbn- 
Then, the possible dominant process would be the SBBN 
process c + d — > a + b if (c, C) and (d, C) are not stable 
against scattering off the background photons. In addi- 
tion, (c, C)+d — > (a, C)+b can be also another dominant 
process if (d, C) is not stable in the thermal bath, for sim- 
plicity assuming (a, C) has larger binding energy than 
(b, C). 11 In these processes, we may expect a Boltz- 
mann suppression factor in the reaction rate e _< 2 CBBN ' /T , 
not e _<3sBBN/T in a similar fashion in SBBN. 

If the SBBN strong interaction a + b — > c + d occurs at 
a shorter time scale than the typical time scale of electro- 
magnetic (EM) interactions of the bound states, we may 
expect that such a short-distance reaction rate should 
not be deviated from the SBBN rate. For D, T, He, Li, 
Be etc, this condition can be realized easily. 

Flux 



the flux has to be estimated by the relative velocity be- 
tween the bound state and the free light element, which is 
controlled by the relative velocity. But even in such cases, 
Vjm ~ Vreac due to the following consideration. Tak- 
ing Vreac = V\ , while the free element goes through the 
target volume, the bound light element rotates with the 
speed V\ = \j2E\>\ n /mx- Then the number of rotations 
would be - V 1 At/2Trr B ~ 0(V 1 /V 2 ) where At - 2r B /V 2 
is the time for the free light clement to go through the 
bound light element, V 2 is the velocity of the free light 
element, and r B is the radius of the bound state. 13 
Then, for the nuclear reaction due to pion exchange, if 
we take V Icac = V\, the flux is the relative velocity V 2 
times 0{V\/V 2 ) which would be ~ V\ 

Next, we consider collisions between a neutral bound 
state a and a neutral (or charged) bound state b. In this 
case, since the target is not a freely propagating parti- 
cle, the speed which controls the flux is not the bound 
light element's V\ <~ V a + V& but the relative velocity V 2 
between the bound states. V 2 is order of the thermal ve- 
locity of the bound state, which is smaller than V\. Then 
V 2 /V\ <~ O(0. 1) at around T=lkeV where neutral bound 
states can be formed. However while the bound states 
collide with each other, the bound element would rotate 
around a CHAMP ~ {V\/V 2 ) times. Therefore even in 
this case, we could estimate y nux <~ V\. 

These considerations imply that we can simply assume 
that the CHAMP in the bound states is a spectator and 



In general, the velocity Vq ux which controls the flux 
might be different from the velocity Keac which controls 
the short distance nuclear reaction. The (crf us io n Vflux) 
would be given by 

(ffusionVflux) = (ffSTT^ + Cp(Vfl ux V reac )....) (35) 

Here, we assume that for the short distance reactions, the 
coefficients, as, up... in CBBN are the same as in SBBN 
12 . Using this approximation, we evaluate the flux. 

First, we consider collisions between a bound state and 
a free light element. Then, once we focus on the 2^2 
collision between the bound and the free light element, 
the relative velocity V\ may be dominated by the speed of 
the bound light element. If we assume that the free light 
element is distributed uniformly in the thermal bath, the 
flux is controlled by V\. On the other hand, in the case 
that the radius of the bound state is smaller than the im- 
pact parameter of nuclear reactions [which is 0(l/m v )], 



11 (c,C)+d — > a+b+C is also possible if it is kincmatically allowed. 

12 If the phase space is modified by the release of a CHAMP after 
the reaction, the difference from the SBBN case would be also 
small if the Q-value is large. 



Corrections for BBN nuclear reaction rates with no threshold 

Here we consider the nuclear reactions containing light 
elements captured by CHAMPs. In this case, as we 
mentioned before, the crucial differences from SBBN are 
in the Coulomb factor and the Boltzmann suppression. 
Since the radius of bound state is very small 0(l/m v ), 
a simple replacement Zx — > Zrx,c) = %x — Zq in the 
Coulomb factor would be a good approximation. How- 
ever the short distance part is also changed because the 
light element captured by a CHAMP has the kinetic en- 



This discussion rely on an assumption that the factorization of 
Coulomb factor and short-distance nuclear fusion is valid. That 
is, we assumed that in the collision, the bound state is not de- 
stroyed before the collision. This would be valid if re ~ If 
the bound state is unstable against incident nucleus, the effective 
^flux/^reac may become smaller than unity. 

Our consideration is based on our approximation that the short- 
distance reaction is the same as that of the SBBN 2 — > 2 process 
between light elements. In the case that the De-Broglie wave- 
length of an incoming nucleus is longer than the Bohr radius 
of the bound state, we may have to solve quantum mechanical 
many-body problems including a bound CHAMP to obtain a 
more reliable result. 
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ergy E^ n not 0(T). As we mentioned before, we as- Then we find 



sume that the short-distance cross section crf nsion takes 
the same functional form of the CM energy as those of 
SBBN, 

Cfusion- Thus the CM energy of the short-distance 
nuclear reaction may be 0{Max{E^ la , E ))- We intro- 
duce these two changes in the estimation of nuclear re- 
action rates. That is, 



X V = (o%+a c P V 2 



)F, 



(aC)b 



{(3) 



2TrZ {aC) Z b a _ 
<r v(V) e 



(36) 



where /3 is the relative velocity between the bound state 
and the incident thermal light clement, and V is the rel- 
ative velocity between the bound light element and the 
incident thermal light element with E — E$. [3 controls 
the amount of penetration in the Coulomb potential. V 
appears in the flux and the short distance cross section. 

Since the short distance cross section would be gov- 
erned by the kinetic energy of the bound light element 
which does not depend on the condition of the thermal 
bath much, the thermal average should be taken only 
for the Coulomb part which implies the evaluation of the 
wave function for an incident thermal light element at the 
position of a bound state. Then the thermal average may 
be taken for the thermal light elements and the thermal 
bound state because the incident thermal light element 
approaches inside the Coulomb field of the bound state, 
not that of bound light elements. We assume that the 
short distance reaction is faster than the EM interaction 
of the bound state. The thermal-averaged cross section 
is calculated as follows. 



(^fusion V ) 

r,2 



9 f 3 

(2ir) 6 n b n {aC) J 

x e 3 e T 



(V)- 



2irZ b Z( aC - ) a 



(3 
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d 3 p r a v(V) 



2irZ b Z( aC * ) a[i( a c)b 



(27r) 3 n r J p r 

2jZ t. Z (aC)°A'(aC)t) Er 

x e pr e T 



8ngTn {aC)b 
(27r) 3 n r 



dyS(yT) Ncw e- {y+ " 



(37) 



where S(y x T) Nev! = a v(V)TrZ b Z ( 
27T2z ! Z (aC)^(aC)ba 2 /T, y = (j/g/4) 1/3 and p (aC)b = 
m (aC) m b/ i m (aC) + TO fc) — m b- Notice that we are assum- 
ing that the short-distance nuclear cross sections have 
the same functional forms of the CM energy as those of 
SBBN. 

Here in the case of Zi aC ^ ^ 0, we relate the new S- 
factor above to the SBBN S factor which could be mea- 
sured by experiments, 

Z {aC) M(aC)6 



S(y x T) Ncw = S(y x T)- 



Z a Pab 



(38) 



(ag sioa V)(T) = 



I 32 E'J 3 S(y x T% 

4V3 ^ {aC )b T2/3 



\ \ T > 



G.U/3 



(39) 



where y x ~ {((» a b/m a )E hin ) + E 9 )/T ~ (E hin + E )/T), 
E G = 2TT 2 Z b 2 Z 2 aC) p (aC)b a 2 and E = Ty . 

For nuclear reaction rates with neutrons like 
7 Be(n,p) 7 Li, since there is no Coulomb suppression or 
Boltzmann suppression if there is no threshold in the 
process [(aC) + n — > (cC) +d], we replace CM energy by 
E C m = ^ab/p(aC)b)E hi n + 3T/2 - E hin + 3T/2 in the 
cross sections because of the change of the kinematics of 
the bound light elements. In addition if the bound state 
is neutral (^( n c) = 0)> * ne Coulomb factor may disappear 
if the bound state is not destroyed before the collision. 
Then the treatments may be similar to the neutron case 
above. Such neutral bound states will be formed in case 
of Z c = 1 (proton, D, and T). 

In the above discussions, we have taken the approxima- 
tions that the light clement is pointlike and does not have 
internal structure, and the selection rules in the nuclear 
reactions are not changed by the trapped CHAMP. 15 



Corrections for BBN nuclear reaction rates with threshold 

Let us consider a case that the SBBN reverse process 
c + d —f a + b has a threshold and the SBBN cross section 
ofa + 6^ c + d can be measured by collider experiments. 



First, 

QsBBN;o 



assuming the condition E\ 



bin(aC):lS 



< 



b.cd 



< E hin(aC y AS + E hin . 1E is satisfied where 



^bin(oC):iS! -E-bin:i£ are the binding energies of the IS 
state of the (a, C) system and of the first excited level 
of the (c, C) system, we can estimate the cross section 
of (c, C) + d — > (a, C) + b by using the information 
of the SBBN forward process a + b — ► c + d. Under 
the above conditions, we may use the detailed balance 
relation on (a, C) + b — > (c, C) + d in a similar fashion 
to the previous discussion. The thermal-averaged cross 
section of (c, C) + d — > (a,C) + b may be written as 
follows. Applying the detailed balance relation and 
the modifications for the forward process which was 



15 We have also assumed the hierarchy between the SBBN strong 
reactions and the EM interactions of the bound states. The 
Bohr radius of CHAMP-light clement system and the typical 
pion-exchange radius would be the same order of magnitude in 
our case, but we can still expect a hierarchy in the coupling 
strengths between EM and strong interaction, which may still 
allow us to factorize the short distance nuclear reaction from the 
effects caused by binding a CHAMP.But if the incoming nucleus 
is very slow, this factorization may break down due to the long 
range nature of the EM force 
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previously discussed, 

( a tusion,(cC)dV) 
_ 9(aC)9b , m b 3/2 

9(cC)9d m d 



'usion, (aC)b 



V)e 



9CBBNaS,cii 



(40) 



where QcBBNab,cd = Q(aC)b,(cC)d- If <3cBBNa6,cd IS Small, 

the Boltzmann suppression might disappear even though 
SBBN has a large Boltzmann suppression. 

For Q C BBN:ab,cd > E hin ( cCy2E where E hin{cC) . 2E is the 
binding energy of the second excited level of the bound 
state, we would not be able to simply apply the detailed 
balance relation for the forward process. But in any case, 
since the crucial point for the processes with a threshold 
is the Boltzmann suppression which comes from the re- 
quirement that the kinetic energy of the incident particles 
overcomes the threshold, if the Q value is smaller than 
that of SBBN, we may expect the milder Boltzmann sup- 
pression in the process, compared to that of SBBN. 16 

In the Zc = 1 case, at a relevant time when capture 
become efficient, the Boltzmann suppression is huge if the 
Q value is O(MeV), and then most of the BBN processes 
are completely decoupled. Hence we ignored the change 
of nuclear reaction rates for the SBBN reverse processes if 
Ocbbn is ~ C(l)MeV, which is a reasonable assumption. 

Next, we consider the reverse process in CBBN, which 
corresponds to SBBN a + b — > c + 7, i.e., (c, C) + 7 — > 
(a, C) + b assuming that the binding energy of (a, C) is 
smaller than that of (b, C) . It is well-known that the 
reaction rate of this forward process is small. Notice 
that the incident photon with the threshold energy of 
the process does not have Coulomb suppression. Thus 
the main origin of suppression is the low abundance of 
the higher energy components of thermal photons. 



(erf • V) 
\ fusion, C7 / 



87T<7 7 

(27r) 3 n r 

(27r) 3 n 7 j , i liliN 



dE r ^ c c) 1 E r( ^ cC - )ry a 0:C7 v(V)e t 



dp 7 p^a 0tC7 v(V)e t 



- ±(^^)3/ 20((afusionXaC)bu))e -^, (41) 



where E. 



p 7 . Although the Q value for the pro- 
cess 3 He(a,7) 7 Be might be smaller than lMcV, the pro- 
cess is negligible at the capture time of 7 Be and this 
reverse process does not seem to provide a significant 
change from SBBN. The change on the threshold en- 
ergy of these photodissociation processes might be im- 
portant when we consider the late-decay effects that the 



16 Here we naively assumed QcBBN:a6.cd > Eo where Eq is the 
Gamow peak energy of the reverse processes. If this condition is 
not satisfied, we may need more careful treatments. 



injected high energy EM energy is thermalized and pro- 
duce a huge number of soft photons, which may destroy 
primordial light elements. 



BBN WITH LONG-LIVED CHAMPS 

Recently WMAP has reported the updated values 
of cosmological parameters under the standard ACDM 
models. We can now check the internal consistency of 
SBBN in the light of WMAP3. It has been pointed out 
that the predicted 7 Li abundance seems too high to agree 
with observed abundances. Also for 6 Li, we have to ex- 
pect an additional production after the BBN epoch, like 
cosmic-ray nucleosynthesis. These tensions or discrepan- 
cies may be tantalizing clues to find new physics. Under 
these circumstances, it is interesting to study the effects 
of new physics. 

In previous sections, we considered the possible 
changes of nuclear reaction rates due to long-lived 
CHAMPs. Here we consider the application for the BBN 
in case of Zc = 1. 



Charged Massive Particle BBN (CBBN) 

We consider the thermal freeze-out of light element 
abundances in CBBN and here we simply ignore the ef- 
fects of possible high-energy injections due to the late 
decay of CHAMPs, which may provide the initial condi- 
tion to consider such late decay phenomenon if the decay 
occurs long enough after the decoupling of the BBN pro- 
cesses. We will later discuss the case where the decays 
occur before the freeze-out. In our estimation, we also as- 
sume the instantaneous captures for each light elements 
at T c = Sbin/40. 17 

In SBBN, abundances of all light elements are com- 
pletely frozen until T ~ 30keV. Since T c is 24kcV for 7 Li 
and 38keV for 7 Be, which is almost the end of SBBN, 
the formations of bound states may change their abun- 
dances. For elements lighter than 6 Li, since the efficient 
captures occur only at below 10 keV, we found that the 
change of nuclear reactions can not recover the processes 
at such a low temperature. This conclusion will hold if 
the difference from our estimation of (crf usion V") is not 
large. Also in most of the reverse processes, the Boltz- 



17 As we mentioned before, if the number density of CHAMPs is 
low tiCHAMp/ n i <S 10~ n , the recombination rate might not 
be sufficiently large compared with the expansion rate of the 
Universe, and we may expect poor captures of CHAMPs. Then 
the most of CHAMPs and light elements will be left as freely- 
propagating ionized particles. Because CHAMPs are supposed 
to decay soon, in this case we can apply the known results in 
decaying particle scenarios in literature. 
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maim suppressions are huge at that time, even though 
we use the new Q value Qcbbn- They do not provide 
any significant change from SBBN. 

Under these circumstances, if the CHAMPs decay be- 
fore the captures of Zx = 1 nuclei, we may expect that 
the sizable change due to the captures occur in elements 
heavier than 7 Li. On the other hand, once the capture of 
proton, D and T starts, since the bound states are neu- 
tral and have no Coulomb suppressions in the nuclear 
reactions, the BBN processes may not freeze out. In the 
next subsection, first of all, we consider the case that 
CHAMPs decay before the captures of Zx = 1 elements 
such as proton, D, or T, which start at below T < 1-2 keV 
(t > 10 6 sec). Later we consider the possible effects due 
to their captures. 

Since the abundances of the light elements differ by 
orders of magnitude, often we can identify the relevant 
processes and neglect the others. For example, when we 
are considering a process a(b, c)d, if n a is much smaller 
than the others (rib, n c and rid), this process is negligi- 
ble for the evolutions of rib, ri c and rid, but important 
only for n a . Therefore elements heavier than 7 Be do not 
significantly affect lighter elements abundances. 



CBBN with Z x > 1 

Here we consider the CBBN with captures of Zx > 1 
nuclei (with Zc = 1). This case will be realized if the 
CHAMP lifetime is shorter than ~ 10 6 sec. In Fig. El we 
show a plot of the light element abundances as a function 
of rj 1 including the corrections only in processes among 
charged light elements (Case A). We can find that the 7 Li 
abundance could decrease much from the SBBN value 
for r/. The decrease is induced by the enhancement of 
7 Li(p,a) 4 He reaction rate due to the capture of 7 Li by 
CHAMPs. As we can see in Fig. [4] the CBBN reac- 
tion rate of 7 Li(p,a) 4 He slowly decreases as a function of 
the energy, compared to that of SBBN at the tempera- 
ture where the Coulomb suppression becomes important, 
which results in later-time decoupling of the process than 
in SBBN. 

We also added processes 7 Be(n, p) 7 Li and 
7 Be(n,a) 4 He, which are associated with neutron 
capture (Case B). In these type of processes, the impor- 
tant change from SBBN is the kinetic energy to be used 
in the nuclear reaction. In SBBN, the typical energy 
is ~ 3T/2. However in CBBN, the energy could be 
0(Ebi n ). If the s-wave partial wave mode dominates the 
process, then the difference might be small. However, 
if higher partial modes such as p-wave dominate, we 
expect significant enhancements of the processes. 

In fact, we found that the change in 7 Be+ 7 Li by the 
modification of the process 7 Be(n,p) 7 Li is negligible. On 
the other hand, since 7 Be(n,a) 4 He is a p-wave domi- 
nant process [87| , the modification of this process should 



be important to predict the primordial abundance of 
7 Be+ 7 Li in CBBN. Unfortunately, we currently only have 
poor experimental data sets for 7 Be(n,a) 4 He. However, 
since there is experimental data for the reverse process 
4 He(a,n) 7 Be [88j, we might be able to theoretically infer 
the cross section of the forward process of 7 Be(n,a) 4 He 
approximately by using detailed balance relations. For 
the moment, however, the experimental data do not have 
sufficient resolutions in the relevant energy region be- 
cause of the significant Coulomb suppression and the 
threshold suppression, to correctly calculate the forward 
rate. Therefore, according to Serpico et al. [52|], as a con- 
servative error we also take a factor of 10 on the process in 
this paper, which does not change the SBBN predictions 
at all and is still consistent with available experimental 
data of the reverse rate [HJ . 

In Fig. 5, we plot the theoretical prediction of 7 Li/H 
(upper panel) and 6 Li/ 7 Li (lower panel) as a function 
of r/. The SBBN predictions are marked by the green 
bands. The red (blue) band is for Case B-I (Case B- 
II) in CBBN. Here we assumed rc c /n 7 = 3.0 x 10~ n 
and the instantaneous capture of CHAMPs. Case B-I 
means that E C m = (fJ>ab/ fi>(aC)b)Ebm + E in a process 
(a, C) + b — > (c,C)+d where we take E to be the Gamow 
peak energy for collisions between two charged elements, 
and to be 3T/2 for collisions between a nucleus and a neu- 
tron. Case B-II means that we take Eqm = -E^in + Eq as 
the CM energy of processes and a 10 times larger value 
of the p-wave part of the cross section of 7 Be(n,a) 4 He 
than that in the standard BBN code [HHB]. In Fig.5, 
it is showed that the modification by a factor of 10 on 
the p-wave partial cross section of 7 Be(n,a) 4 He does not 
change the SBBN prediction (Case B-I) but must be im- 
portant in CBBN (Case B-II). 

We have also checked the reverse process of 
10 B(p,a) 7 Be. The threshold in this process can become 
smaller, which may induce milder Boltzmann suppression 
than that of SBBN. However, we found that this rate is 
simultaneously suppressed strongly by the Coulomb fac- 
tor, and therefore this effect is irrelevant. 

Finally we warn the readers again that our results 
rely on the assumption that the short-distance nuclear- 
reaction rates have the same functional form of the CM 
energy as those of SBBN. In addition, we assume that 
by relevant elements, the energy to excite nucleons into 
higher levels and the binding energy by a CHAMP are of 
the same order of magnitude. To obtain a quantitative 
conclusion, further efforts to estimate the errors in the 
short-distance nuclear reaction rates must be important. 
For example, in 7 Be(n,p) 7 Li, the change of the nuclear 
reaction rate can directly affect on the final abundance 
of 7 Li(= 7 Li+ 7 Be). However notice that well before the 
elements lighter than Li are captured by CHAMPs, the 
SBBN processes are completely decoupled. Even though 
the errors induce larger reaction rate, if it were within 
an order of magnitude level, nuclear reaction would not 
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FIG. 3: Theoretical predictions of Y p , D/H, 3 He/H, 7 Li/H 
6 Li/H, 3 He/D and 6 Li/ 7 L i as a function of r) in standard 
BBN (green) and CHAMP BBN in case A (red). Here we have 
assumed the instantaneous capture of CHAMPs and nc /n 7 = 
3.0 x KT 11 . 



overcome the expansion rate again, and our conclusion 
would not be changed, because the Coulomb suppression 
is significant, and the neutron abundance is very small. 



FIG. 4: Ratio of nuclear-reaction rates of SBBN and CBBN 
in Case B-I and B-II as a function of the cosmic temper- 
ature for the relevant processes. Here we assumed the in- 
stantaneous capture of CHAMPs by the nuclei. Case B- 
I means that Ecm = (flab/ 'fJ>( a c)b)Ebin + E in & process 
(a, C) + b — > (c, C) + d where we take Eq to be the Gamow's 
peak energy for collisions between two charged elements, and 
to be 3T/2 for collisions between a nucleus and a neutron. 
The Case B-II means that we take Ecm = -Ebm + Eo as the 
CM energy of processes and 10 times larger value of the p- 
wave part of the cross section of 7 Be(n,a) 4 He than that in the 
standard BBN code El, El. 



CBBN with Z x = 1 

Next we discuss the possible effects due to captures 
of Zx = 1 nuclei. Since the bound states are neu- 
tral, the nuclear reactions in TABLE II may not have 
Coulomb suppression and might be significantly changed 
from those of SBBN. 18 

We consider the case that Zc = 1 nuclei are 
captured instantaneously at temperatures below each 
T c (~ 0(1) keV). The captures of T only provide a 
significant change in T itself and 7 Li even if we as- 
sume the instantaneous captures because of their poor 
abundances. The captures of D result in large enhance- 
ments for the processes listed in TABLE II. In partic- 
ular, since T(d,n) 4 He, 3 He(d,p) 4 He, 7 Li(d,na) 4 He and 
7 Be(d,pa) 4 He have large cross sections, the reaction rates 
may be able to become larger than the expansion rate 
again at a later time. Their decoupling does not occur 
soon because of the absence of the Coulomb suppres- 



Notice that here we simply assumed that the bound state is not 
significantly disturbed before the nuclear fusion reactions. 



sions. If the captured D abundance is larger than 3 He, 
D and 3 He mainly burn into 4 He through 3 He(d,p) 4 Hc. 
Then the abundance of 3 He can decrease, and the abun- 
dance of D becomes close to rip — ns He . On the other 
hand, however, D(d,p)T and D(d,n) 3 He do not change 
the abundance of D so much. In addition, T(d,n) 4 He, 
7 Li(d,na) 4 He, and 7 Be(d,pa) 4 He do not a change the 
abundance of D either, but might decrease the abun- 
dances of T, 7 Li and 7 Be because of the small abundances 
compared with that of D. 

Although the process D(a,7) 6 Li has very small re- 
action rate in SBBN since the abundances of the inci- 
dent particles (D and 4 He) are sufficiently large, this 
process can produce large amount of 6 Li. 19 The cap- 
ture of protons can reduce the 6 Li and 7 Li abundances 



Recently Pospelov pointed out that the cross section of 
D(o,7) 6 Li might be significantly enhanced by considering the 
virtual photon absorption due to a bound CHAMP 89], which 
might significantly overproduce 6 Li. Since his paper appeared af- 
ter the completion of this work, we have not included this effect 
in this paper. 
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FIG. 5: Theoretical prediction of 7 Li/H (upper panel) and 
6 Li/ 7 Li (lower panel) as a function of the baryon-to- photon 
ratio. The SBBN predictions are marked by the green bands. 
The red (blue) band is for Case B-I (Case B-II) in CBBN. 
Here we assumed nc/fiy = 3.0 x 10 -11 and the instantaneous 
capture of CHAMPs. The definition of Case B-I and Case 
B-II are same as those in Fig. [4] 



through 7 Li(p,a) 4 He and 6 Li(p,a) 3 He. On the other 
hand, there are no significant changes on D, T and 7 Be 
abundances because the associated processes are radia- 
tive ones, which are relatively suppressed. 

In the relevant epoch for the captures of p, D and T 
(T < 0(1) keV), the condition to overcome the expan- 
sion rate is that the reaction rates are larger than that of 
10 4 cm 3 /sec multiplied by the captured number density 
of Up. Then in processes T(d,n) 4 He, 7 Li(d,na) 4 He and 
7 Be(d,pa) 4 He, even if the decrease of reaction rates were 
within a factor of O(10) due to some ambiguities such as 
capture rate of D, we could still expect the decrease on 
7 Li and 7 Be abundances. As we showed before, since the 
changes on light element abundances by the captures of 
Zx > 1 nuclei might be small, the initial condition of 
light element abundances for such a later-time CBBN by 
captured Zx — 1 nuclei might be the same as those of 
SBBN. However notice that the above conclusions rely on 
the number density of the captured Zx = 1 nuclei very 
much. If the number density of CHAMP is not large, 
the captures weaken, and the changes become milder. 
For example, taking possible capture fractions, O(10 -5 ), 
O(0. 1), and O(10~ 2 ) for proton, D, and T, respectively, 
we show the results in Fig. [6j In such cases, the nuclear- 
reaction rates for 7 Li, 6 Li and 7 Be become more rapid 
than the expansion of the universe, and we expect that 
7 Li and 7 Be decrease without changing D, 3 He and 4 He 




FIG. 6: Theoretical predictions of 7 Li/H and 6 Li/ 7 Li as a 
function of r\ in SBBN and CBBN with their theoretical er- 
rors at 95 % C.L. Here we took the fractions of the captured 
proton, D, and T to be 10 -5 , 10 _1 , and 10 -2 , respectively. 
The WMAP value of r\ at 95 % C.L. is also indicated as a 
vertical band. We can find that the primordial values of 6 Li 
and 7 Li in the CBBN may be in the range of the observed 
abundances, which may simultanously solve current 6 Li and 
7 Li problems pointed out in SBBN. 



abundances. The 7 Li abundance is determined by the 
competition between two processes, 7 Li(p,a) 4 He for pro- 
ton capture and T(a,7) 7 Li for T capture. The 6 Li is 
controlled by the production reaction D(a,7) 6 Li, and the 
destruction reaction 6 Li(p,a) 3 He. In the case of Fig. [51 a 
sizable amount of 6 Li is produced, and the predicted pri- 
mordial value of 6 Li/ 7 Li approximately agrees with the 
observational data without assuming any chemical evo- 
lution scenarios. 

On the other hand, notice that some ambiguities might 
still exist in the nuclear-reaction rates. For Zx = 1 nu- 
clei, because the bound state with Zc = 1 CHAMP has 
a larger Bohr radius than those of Zx > 2 nuclei, the 
electromagnetic disturbance on the bound state before 
the nuclear fusion reactions occur would have to be more 
carefully considered. If the bound state is electromag- 
netically destroyed by an incident heavier nucleus, the 
factorization of Coulomb part and short-distance nuclear 
reaction part does not work well, and the nuclear-reaction 
rate may be changed from the value of our calculations. 20 



For the collision of neutral bound state with Zx = 1 nucleus, the 
formation of molecule may be important to evaluate the nuclear- 
reaction rate. 
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process (p,C) 


The reaction rate (cm /sec/mol) 


7 Be((p,C), 7 )( 8 B,C) 


(4 - 24) x 10 2 


7 Li((p,C),a)( 4 He+C) 


3 x 10 6 


6 Li((p,C),a)( 3 He+C) 


1 x 10 8 


6 Li((p,C), 7 )( 7 Be+C) 


3 x 10 3 


T((p,C), 7 )( 4 He+C) 


2 x 10 2 


D((p,C), 7 )( 3 He+C) 


40 


process(D,C) 


The reaction rate (cm 3 /sec/mol) 


(D,C)(p, 7 )( 3 He+C) 


40 


(D,C)(Q, 7 )( 6 Li+C) 


0.6 


(D,C)(d,n)( 3 He+C) 


7 x 10 6 


(D,C)(d,p)(T+C) 


4 x 10 6 


T((d,C),n)( 4 He+C) 


1 x 10 9 


3 He((d,C),p)( 4 He+C) 


4 x 10 8 


7 Li((d,C),na)( 4 He+C) 


1 x 10 8 


7 Be((d,C),pa)( 4 He+C) 


3 x 10 8 


processes(T,C) 


The reaction rate (cm 3 /sec/mol) 


(T,C)(p,7)( 3 He+C) 


2 x 10 2 


(T,C)(d,n)( 4 He+C) 


1 x 10 9 


(T,C)(a, 7 )( 7 Li+C) 


2 x 10 3 



TABLE II: Table of typical values for nuclear-reaction rates 
with a captured proton, D, and T. Here we ignored nuclear- 
reaction rates related to the neutron because they are negligi- 
bly smaller for T < T c . Here we also ignored processes which 
include elements heavier than 7 Be in the initial state. 



Also if a large amount of CHAMPs survive until such 
late times (> 10 6 sec), we would have to simultaneously 
consider both effects due to the captures by the Zx = 1 
nuclei, and subsequently the EM energy injections by the 
decaying CHAMPs at a later time. 



Late decays of long-lived CHAMP 

We have discussed the change of light element abun- 
dances before the decay of CHAMPs. On the other hand, 
the decaying CHAMPs might induce additional changes 
of primordial light element abundances, which have been 
studied by several groups. The effects highly depend on 
the decay products, i.e., electromagnetic or hadronic cas- 
cades [45l 

10 4 sec, 



I46J, |4& |49| 

At an earlier epoch before t = 10 4 sec, only the 
hadronic energy injection is important, and there is al- 
most no constraint from EM energy injections. There- 
fore, at such an epoch, even though the injected energy is 
not small, if the branching ratio into hadronic cascade is 
sufficiently suppressed, there are no significant effects on 
the primordial light element abundances. Such a case is 
well known if CHAMPs decay into leptons with a branch- 
ing ratio into hadrons of the order of 0(1O~ 3 )-0(1O -6 ). 
For late decays after 10 4 sec, the amount of energy release 
may be highly constrained by the EM energy injections. 



In the following sections, we consider possible new 
changes by taking into account the capture of CHAMPs. 



Are there corrections on the evaluation for the primary 
energy injection by CHAMP decays ? 

Since the binding energy of CHAMP bound states is 
below the nuclear binding energy of the light elements, 
the recoil of nucleons inside the captured light element 
due to CHAMP decays would not destroy the light el- 
ement. On the other hand, the decay products of the 
bound CHAMP might directly hit the bound light ele- 
ment and destroy it. Let us consider the case that the 
primary decay product is a charged lepton as an exam- 
ple. Of course, if the lepton is a tau, the tau lepton soon 
decays into hadronic particles. However, the lifetime of a 
tau lepton is long enough to go through the Bohr radius 
of the bound state. Thus we will deal with all kinds of 
leptons in a similar fashion. 

Naively we may speculate that the light elements are 
distributed inside the radius r x x ~ A^l^m^ and the 
CHAMP stays somewhere inside the radius. Then the 
number density of quarks (or nucleons/ nuclei) inside a 
bound nucleus is roughly 



"boundx ~ Mr— f - m%. 
rx 

The mean free path is roughly estimated by 

1 



A 



mfp 



O X TlboundX 



(42) 



(43) 



where we have chosen a ~ 2iTa 2 ft where t is the Man- 
delstam variable t for the momentum transfer from a pri- 
mary decay product (a charged lepton) to a bound light 
element. Then, the naive probability of the primary de- 
cay product (a charged lepton) scattering off a quark (or 
nucleon/nucleus) inside the nucleus is 



Prob 



2r x 



A 



mfp 



O(10" 



(10GeV) ; 
t 



A l/3 



(44) 



Among light elements, 4 He destruction would be most 
dangerous. If we assume that all of 4 He's are completely 
captured by CHAMPs, if such a probability is below 
10~ 4 , the change on D/H, 3 He/H abundances due to the 
direct collision will be below the O(10 -5 ) level, which 
may not disagree with observed abundances. Elements 
heavier than the destroyed parent nuclei should not be 
directly produced significantly. 21 



21 However, in a recent work 9(J, they have pointed out the pos- 
sibility that energetic T and 3 He which are produced from the 
destruction of the bound 4 He can nonthermally produce sizable 
amount of 6 Li. 
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We can find that, if the momentum transfer from pri- 
mary decay-product is hard (t > (lOOMeV) 2 ) 22 , the light 
element bound by a CHAMP is sufficiently transparent 
and may not disturb the SBBN prediction for elements 
lighter than 4 He. In this case, for the evaluation of 
the primary energy injection, past studies in the liter- 
ature will be a good approximation, which considered 
that CHAMPs are freely propagating in a thermal bath. 
For the secondary products through the hadronization of 
a recoiled quark or direct production of nucleons/nuclei, 
the above probability will be identified as the hadronic 
branching ratio for a CHAMP decay, which may provide 
only negligible effects on elements lighter than 4 He. On 
the other hand, we can consider another extreme case 
where the momentum transfer is sufficiently soft. For 
example, if the energy is smaller than the nuclear bind- 
ing energy, the charged lepton of decay products could 
not inelastically scatter off the bound light element. In 
the middle range between them, we may have to simul- 
taneously consider the direct collision and EM/Hadronic 
cascade induced by the late decay which was considered 
before. 

In the case of hadronic decays, we may replace a by 
the strong coupling a s in the above estimation. Then, 
we find that for a sufficiently hard momentum transfer 
t > (lGeV) 2 , the bound light element is still transparent. 



Other new possible corrections to BBN with late time energy 
injection 

There are three types of other possible effects on light 
element abundances by the late-time decaying CHAMPs 
when some fraction of such light elements are captured 
by CHAMPs. 

The first type originates from the change on the 
Coulomb barrier and the kinematics of background light 
elements as targets for nonthermal processes by their 
own bound states. This change could be important for 
the hadronic-decay scenario. The injected high-energy 
hadrons eventually lose their energy due to the ther- 
mal interactions and become nonrelativistic and collide 
with the background light elements. If the target nuclei 
are captured by CHAMPs, the reaction rates for various 
hadronic processes might be different from experimental 
values. 

The second type is related to the change of the Q value. 
This might be important for both high-energy hadronic 
and EM energy injections, especially in high-energy pho- 
ton injections. The injected high-energy photons pro- 



duce many soft photons through the EM cascade before 
the scattering off background light elements. Then the 
spectrum of the soft photons has a cutoff at the energy 
above the threshold of electron-positron pair creation, 
which depends on the cosmic time or the cosmic temper- 
ature. Only when the cutoff energy is higher than the 
threshold energy of the photodissociation are the target 
nuclei destroyed. The change in the Q value may mod- 
ify the epoch when the light element is destroyed by the 
photodissociation processes. 23 

The third type is related to neutron injection from late 
decays. If the decay occurs while some of the SBBN pro- 
cesses are still active, the high-energy hadronic injections 
might produce many neutrons. At late time (t > 100 
sec), since neutrons have an extremely low abundance 
due to the (3 decay, the produced neutrons can signif- 
icantly affect the light element abundances because the 
related nuclear reactions do not have a Coulomb suppres- 
sion. The neutron injection at around 10 3 sec was dis- 
cussed as a solution to obtain low 7 Be abundance by the 
destruction of SBBN 7 Be through 7 Be(n,p) 7 Li and sub- 
sequently 7 Li(p,a) 4 He 28|, |4y, |9l| . In our scenario, there 
may exist some differences from the previous studies. As 
we mentioned before, in CBBN, 7 Be(n,a) 4 He could be 
more important for the 7 Be abundance than 7 Be(n,p) 7 Li 
because the center-of-mass energy in the process can be 
completely different from that of SBBN. Thus the neu- 
tron injections from the late decaying CHAMPs may en- 
hance the destruction ability of 7 Be, and the effects could 
be different from the noncaptured case. Notice also that 
there may still exist unknown errors even on the reaction 
rate of 7 Be(n,p) 7 Li with captured 7 Be, as was discussed 
before. 

The modification of the reaction rate of 3 He(d,p) 4 He 
might be interesting with decaying CHAMP scenario be- 
low IkeV. If a sizable fraction of D is captured and the 
reaction rate of the process is enhanced due to the cap- 
ture of a CHAMP, the destruction rate of 3 He might be- 
come more rapid than the production due to late de- 
cays of CHAMPs, which may weaken the bound on 3 He 
production due to the decay. This new possibility may 
relax the 3 He/D bound, which is generally the most se- 
vere constraint on radiatively or hadronically decaying 
massive-particle scenarios at t > 10 7 sec. 

Considering the above possibilities, it is important to 
reanalyze the effect of the late-time decaying CHAMPs 
in their bound states with the light elements 86]. 



The energy transfer due to the momentum transfer < (lOOMeV) 2 
may be below the typical threshold ~ 0(10) MeV to destroy a 
bound light element by NR nucleon/nuclei scattering inside the 
light clement. 



On the other hand, we may also have to take care of the destruc- 
tion of bound state due to huge soft photons from high energy 
photon injection. 
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DISCUSSIONS 

If the CBBN prediction is not be significantly dis- 
turbed by late decays, the superWIMP dark matter sce- 
nario would be interesting. Here we discuss how much 
relic of CHAMPs might be allowed in this scenario. Since 
the number density of CHAMPs is important to evaluate 
the capture rate of light element, we consider the possible 
constraints on the number density of CHAMPs, assum- 
ing that the whole dark matter originates from the two 
body decay of CHAMPs into a dark matter and a SM 
particle. We consider the free streaming by the whole 
dark matter produced from CHAMP decays. The relic 
density of CHAMP is 



n c = 



m c 

771 DM 



n 



DM 



(45) 



As we found before, the capture rate is governed by the 
number density of the CHAMP. 



77C 



77DM 



3 x 10" 



L1 100GcV » DM 
777DM 0.23 



(46) 



Hence the lighter mass of the dark-matter allows larger 
CHAMP abundance. Since keV warm dark matter is still 
allowed from Lya data 92J, we naively require that the 
dark matter is nonrelativistic at T =keV. Then we find 
a following condition. 



u < 1.0 



10 6 sec 



t 



(47) 



where u = y/\p i Pi\/m and pt is the three momentum 
of dark matter. Assuming the two body decay, the 
four velocity at the decay time is u = (fftcHAMP — 
777 dm )/2?77dm777champ- Then we find that, for lifetime 
~ 10 4 sec, u ~ 20 may be allowed. Then it is possible to 
take 77champ/^7 ~ O(10 -9 ), which will lead to a con- 
siderable capture rate. 

For the case that the decaying CHAMPs contribute 
only to part of the dark matter, or their contribution is 
negligible, the above constraint may not be applicable. 



CONCLUSION 

In this paper, we have discussed the role of long-lived 
charged particle during/after the BBN epoch. We found 
that the existence of CHAMP during the BBN epoch can 
change the light clement abundances if the capture rate of 
CHAMP by light elements is sufficiently large. Since the 
bound state for heavier elements tends to be more stable 
against the destruction by the background photon, the 
abundances are modified only for heavier elements such 
as Li and Be, thanks to the capture at an earlier time be- 
fore the nuclear reactions decouple. On the other hand, 
the abundances of lighter elements such as D, T, 3 He, 



and 4 He are unchanged. In fact, even though more work 
needs to be done to find quantitative results, we have 
shown that the capture of CHAMPs may possibly have 
some impact on the BBN prediction of the primordial 
7 Li abundance. Our approach to consider the cosmolog- 
ical effects of the formation of the CHAMP bound states 
should also be attractive in some particle physics mod- 
els [alii. 

To understand CBBN more correctly, we need to un- 
derstand the nuclear fusion rates and the capture rates 
more precisely. However, unfortunately there are still 
some uncertainties in the experimental data of the reac- 
tion rates at present. We expect that the future nuclear 
experiments will clarify these points. If future collider 
experiments find a signal of long-lived charged particle 
inside the detector, the measurement of lifetime and de- 
cay properties of the charged particle will provide new 
insights to understand the phenomena in the early uni- 
verse in turn. 
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